Experimental observation of the 'Tilting Mode' of an array of vortices in a dilute 

Bose-Einstein Condensate 
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We have measured the precession frequency of a vortex lattice in a Bose-Einstein condensate of 
Rb atoms. The observed mode corresponds to a collective motion in which all the vortices in 
the array are tilted by a small angle with respect to the z-axis (the symmetry axis of the trapping 
potential) and synchronously rotate about this axis. This motion corresponds to excitation of a 
Kelvin wave along the core of each vortex and we have verified that it has the handedness expected 
for such helical waves, i.e. precession in the opposite sense to the rotational flow around the vortices. 
The experimental method used to excite this collective mode closely resembles that used to study 
the scissors mode and excitation of the scissors mode for a condensate containing a vortex array 
was used to determine the angular momentum of the system. Indeed, the collective tilting of the 
array that we have observed has previously been referred to as an 'anomalous' scissors mode. 
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Since a vortex was first nucleated in a dilute Bose- 
condensed gas jl| , there has been a considerable effort to 
understand the dynamical behaviour of individual vor- 
tices and vortex arrays. This research has, in part, been 
driven by the parallels between dilute gas systems and 
more complex superfluids such as He^, however with the 
Bose-condensed gases it has proved to be straightforward 
to obtain images of vortices and to measure properties of 
a single vortex. The precession of a single vortex has been 
investigated theoretically y, |3| , and studied experimen- 
tally in a nearly spherical Bose-condensate [j] . It is also 
possible to nucleate many vortices in a Bose-condensed 
gas and these form a regular Abrikosov lattice [5|, |a, Q] ■ 
The work described here can be described in two comple- 
mentary ways: (a) as an extension of the previous work 
on the precession of a single vortex to the case of an ar- 
ray of vortices in an anisotropic trap where the collective 
motion of the vortices is relatively rapid, or (b) as the 
excitation of the lowest-energy Kelvin wave of a vortex 
lattice. This second viewpoint is described in more detail 
below. 

Vortices break the degeneracy of certain modes in the 
normal Bogoliubov excitation spectrum for a trapped 
condensate. This splitting has been observed for both 
the m=2 quadrupole mode '&' and also for the scissors 
mode |2|; in the latter case the precession that arises 
when the condensate has some angular momentum leads 
to a 'superfluid gyroscopic' motion |lO|. More recently an 
excitation of a vortex in the form of a helical Kelvin wave 
has been detected [ll| ; collective oscillations of a vortex 
lattice called Tkachenko modes have been observed and 



their frequency measured 
paper Chevy and Stringari 



. In a recent theoretical 
have extended the hydro- 



dynamic theory of a Bose-condensed gas to describe an 
array of vortices and they predict that an array precesses 
in a similar way to a single vortex, i.e. all the vortices si- 
multaneously undergo a motion equivalent to the lowest- 



energy Kelvin wave of a single vortex (see Fig. ^ . (In 
these waves the orientation of the vortex core changes 
and, unlike the Tkachenko modes, the Kelvin waves do 
not arise in a two-dimensional system.) In this paper 
we present measurements of the frequency of this 'collec- 
tive tilting' mode of the vortex array by an experimental 
method similar to that used to study the scissors mode of 
the condensate [lj|. Indeed Chevy and Stringari []j| re- 
fer to the collective tilting of the array as an 'anomalous' 
scissors mode. Their hydrodynamic theory predicts that 
the frequencies of this mode and the two normal scissors 
modes, in a reference frame rotating at frequency Oq, are 
given by the roots of the cubic equation: 



± 2nQiJ - uj{luI 



-uji-nf,)T'2^o^i = 0, 



(1) 



where ujj_ and lUz are the radial and axial angular oscil- 
lation frequencies, fig = nvh/2m is the effective rotation 
frequency of the condensate; n„ is the number of vor- 
tices per unit area and m is the mass of the condensed 
isotope. The hydrodynamic theory assumes a uniform 
distribution of vortices within the condensate so that it 
mimics the rotation of a rigid body. 

We denote the (angular) frequencies of the upper and 
lower scissors modes as w^ and w/ respectively and the 
frequency of the tilting mode by ujt- These frequencies 
are calculated in the non-rotating, laboratory reference 
frame. Generally we shall express frequencies in terms 
of their fraction of the radial frequency of the trap. The 
three solutions of Eq. ^in a potential with u)z/uj±_ = v8 
are shown in Fig [21 (a) . This shows that when JIq = 
the scissors mode frequency is three times the radial trap 
frequency, and that the splitting between the upper and 
lower scissors modes is about 2r2o — "^^t- Fig. [21(b) shows 
a plot of the three quantities aj„/(jjj^ — 3, 3—u!i/uj±, ujt/uj± 
and also the ratio (a;„ — ijji)/2uif This shows that the 
frequencies of all the modes vary approximately linearly 
with r^o (angular momentum); in particular lo^ is very 



close to linear over the entire range but both the lower 
scissors mode and the tilting mode show a noticeable de- 
viation. The given ratio has a value within 15% of unity 
throughout the range. Note that the frequency of the 
tilting mode tends to zero as Qq tends to zero indicating 
that this is a mode of the vortices themselves, in contrast 
the normal scissors modes are modes of the condensate 
whose frequency is perturbed by the presence of vortices 
(angular momentum). It can be seen directly from the 
cubic equation that for a condensate in a spherically sym- 
metric potential {u>x = ^y = ^z) the 'tilting' mode has 
zero frequency; this arises because vortices have the same 
energy for any orientation in a spherical cloud and so if 
the orientation of a vortex is changed with respect to 
some axis then it will simply remain at the new angle. 

Our experiment uses evaporative cooling in a TOP 
(Time-averaged Orbiting Potential) trap to create a 
Bose-condensate that contains ~ 10^ atoms of the Rb^^ 
isotope in the \F — l,mp — —1) state. The oscillation 
frequencies of atoms in the trap are 61 Hz radially and 
172.5 Hz along the z (vertical) direction. After a conden- 
sate has been formed we rotate the trapping potential 
around the z-eods at a frequency of 46.5 Hz, and adia- 
batically change from a cylindrical trapping potential to 
a rotating elliptical potential over a period of 200 ms. 
The potential is rotated for 700 ms at a final elliptic- 
ity where Wx/wy = 0.95, before it is ramped back to an 
axially-symmetric (circular) potential. During this pro- 
cess vortices nucleate at the edge of the condensate and 
then move towards the centre of the trap [i3| . This proce- 
dure introduces ten or more vortices into the condensate 
(Fig.EJ. 

The trapping potential is tilted by the addition of an 
oscillating field in the z-direction, in phase with the usual 
rotating field in the TOP trap. When no vortices are 
present, tilting the trapping potential excites the scissors 
mode as described in jlj|. This scissors motion is the 
superposition of two degenerate modes characterized by 
functions of the form f{r)xz and f{r)yz [ig- The pres- 
ence of one or more vortices, breaks the degeneracy and 
leads to modes described by f{r)z{x±.iy); these are eigen- 
states of the z-component of orbital angular momentum 
Lz ■ A sudden tilt of the trapping potential excites a su- 
perposition of these two counter-rotating modes and the 
resulting scissors oscillation precesses at a rate equal to 
the freque ncy splitting between the two modes divided 
by two Isl, llll . This motion is illustrated in Fig. ^ for a 
condensate that contains an array of many vortices. In 
this case the effective rotation frequency of the conden- 
sate r^o is much greater than that for a single vortex. A 
fit to these data gives the frequencies of the upper and 
lower scissors modes as /„ = cOu/'^tt = 211.8 ± 2.0 Hz and 
/; = 156.3 ± 2.0 Hz respectively. We find that this corre- 
sponds to an average angular momentum per particle of 
(lz) = 8.4±0.4 n, for a total number of atoms TV = 75 000. 
A vortex at the trap centre contributes h of angular mo- 



mentum per particle but off-centre vortices have a smaller 
contribution '2|. From these measured frequencies and 
Eq. n an effective rotation of Qq = 0.495 ± 0.019 can 
be deduced, which implies a tilting mode frequency ujt = 
23.5 ± 0.9 Hz according to the hydrodynamic theory. 

To excite the tilting mode of the vortex array the tilt 
angle of the trapping potential (relative to the z axis) was 
driven at a frequency of 61 Hz for two complete cycles of 
oscillation with an amplitude of 0.07 radians. After exci- 
tation, the condensate was held in the trap for a variable 
amount of time during which the tilting mode evolved at 
its natural frequency; the condensate was then released 
and allowed to expand freely for 19 ms before a laser 
beam was flashed on to record an image of the cloud. 
(This imaging was destructive.) When no vortices were 
present in the condensate the response to the driving 
was an excitation of the normal scissors mode with a 
small amplitude, as in Fig. El because this mode is far 
from resonance. For a condensate containing vortices, 
however, there was near-resonant driving of the tilting 
mode, which resulted in large amplitude oscillations of 
this mode. The tilting of the vortex array leads to a 
change in angle of the condensate and the projection of 
this motion onto the imaging direction gives a signal sim- 
ilar to that of the normal scissors mode but at a much 
lower frequency. In a random sample of the absorption 
images the vortex cores line up with the imaging beam: 
in these images it is possible to see that the vortices are 
indeed tilting in unison when the mode is excited (see 
Fig. EJ. The tilting mode was excited both along the 
direction of the imaging beam and in a direction perpen- 
dicular to it (Fig. IT)). This allowed us to determine the 
direction of rotation of the tilting mode, and to verify 
that it rotates in the opposite direction to the initial ro- 
tation that creates the vortices (the initial rotation has 
the same direction as the flow around the vortices). In 
a Kelvin wave the vortex core has the form of a helix of 
a particular handedness with respect to the direction of 
rotation of the vortex; this property also applies to the 
collective tilting of the vortex array. A sinusoidal fit to 
the data gives a frequency of 57.7 ± 1.3 Hz for the tilting 
mode. 

The number and position of the vortices varies from 
shot to shot resulting in slightly different frequencies of 
the tilting mode so that after a few cycles the observa- 
tions have more fluctuations (i.e. this method of record- 
ing the data has a dephasing analogous to the transverse 
relaxation in magnetic resonance techniques, and we have 
not measured the damping of this mode equivalent to lon- 
gitudinal relaxation time.) The uncertainty produced by 
variations in the initial conditions could be reduced by 
taking more data for each evolution time and averaging, 
as in previous work on the superfluid gyroscope [3|- 

Using the measured values of /„, /; and ft we find 
the ratio of frequencies (/„ — fi)/{2ft) = 27.7/57.8 = 
0.48 ± 0.03 that is not consistent with predictions of the 



hydrodynamic theory, shown in Fig. |5] (b). If we as- 
sume that the scissors modes have the correct frequen- 
cies, then the measured value for the tilting frequency is 
about 2.5 times the hydrodynamic prediction. We should 
emphasize, however, that we cannot determine from the 
frequency measurements whether the scissors modes are 
more accurately described by the theory than the tilting 
mode. One way of checking this is to estimate the con- 
tributions to the angular momentum from the number 
of observed vortices and their position within the cloud, 
and compare this with the value obtained from the scis- 
sors mode splitting: a vortex at radius r in a condensate 
of Thomas-Fermi radius R contributes Nh{l — r^/i?^)^/^ 
to the total angular momentum. In this way we found a 
total angular momentum per particle of 8.5 ± 1.0 fi for 
the image in Fig. 13 which is in good agreement with the 
value calculated from the scissors mode frequency. 

If the prediction of the scissors mode frequency is cor- 
rect, then we need to explain the difference between the 
predicted and measured values of the tilting mode fre- 
quency. Firstly, to eliminate the possibility that we were 
exciting a higher order Kelvin mode, the tilt angle of 
the trapping potential was driven at 25 Hz after a vor- 
tex array had been nucleated. This driving frequency is 
close to the hydrodynamic prediction, but no response 
from the condensate was observed. If we are definitely 
driving the lowest order mode, the case for which is sup- 
ported by the straightness of the vortex cores in Fig. [SJ 
then the difference is most likely explained by considering 
the range of applicability of the hydrodynamic theory; 
indeed, the theory relies on the fact that it is possible 
to average physical quantities over domains containing 
several vortices |l3j . For a single vortex a full numeri- 
cal simulation is currently possible, as shown in Fig. ^ 
for our experimental conditions. The tilting mode fre- 
quency is calculated to be ft — 0.27ujz, where tUz is the 
axial trapping frequency; the splitting between the scis- 
sors modes is calculated to be 0.004^^, thus for a single 
vortex (fu — fi)/{2ft) — 0.09. This deviates even further 
from the hydrodynamic prediction that the frequency ra- 
tio is about unity than the results for our small vortex 
array. The hydrodynamic theory only gives an accurate 
description for arrays that contain a large number of vor- 
tices, the vortex arrays in this experiment do not contain 
a sufficient number of vortices to be in the regime where it 
is possible to average over domains of vortices and apply 
the hydrodynamic theory (Fig. QJ ; increasing the value 
of experiment in this intermediate region. The hydrody- 
namic theory should, however, be a better approximation 



for larger vortex arrays as in references |a, |7| ■ 

In conclusion, this experimental measurement of the 
frequency and direction of rotation of the collective tilt- 
ing mode of the vortex array qualitatively supports the 
predictions made by Chevy and Stringari JI3 . The cause 
of the discrepancy between our measurements and the 
hydrodynamic predictions could be investigated exper- 
imentally with a larger condensate that can contain a 
higher number of vortices or theoretically by a numerical 
simulation as in |17|. The experimental method for di- 
rect excitation of a vortex (or vortex array) can also be 
used to study higher-order Kelvin waves. 
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FIG. 1: The lowest-energy Kelvin wave of a single vortex 
corresponds to a rotation of the tilted vortex line in the 
opposite direction to the rotation of the condensate. The 
plots show one complete rotation of the vortex going from 
top to bottom. In the tilting mode of a vortex array ex- 
cited in this experiment each vortex synchronously undergoes 
a similar motion. The plots show surfaces of constant den- 
sity |<I'p=9e-4 of a condensate cloud (wave function normal- 
ized to one) from a numerical simulation using the Gross- 
Pitaevskii equation. The initial state is a single centred vortex 
(aligned along the vertical axis); this is excited by resonant 
driving for two vortex precession periods by tilting the con- 
fining potential about an axis in the horizontal plane, as in 
the experiments, with an amplitude of 4°. In this simulation 
iV = 75000, f± = 61 Hz, /^ = 172.5 Hz. 
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FIG. 2: Calculation of the mode frequencies predicted by Eq. 
Q for a trapping potential with uiz/ijj±_ — \/8. (a) uju/loi. 
(dot-dash line), uoi/uj± (dotted line), ujt/u}± (solid line), (b) 
oju/ojx — 3 (dot-dash line), 3 — loi/uj±^ (dotted line), uit/u}i_ 
(solid line). The plot also shows the ratio (uju — oJi^jliOt 
(dashed line). 




FIG. 3: Absorption image of a vortex array, taken along z 
axis. 
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FIG. 4: Exciting the normal scissors mode in the presence of 
vortices produces a precession that can be used to infer the 
angular momentum of the vortex array. 
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FIG. 5: (black squares) Off resonant excitation of the scissors 
mode when no vortices are present, (white circles) Resonant 
driving of the 'tilting mode' in the presence of a vortex array. 




FIG. 6: Absorption images of the vortex lattice, viewed in a 
radial direction. 
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FIG. 7: (black squares) Anomalous mode initially excited per- 
pendicular to probe beam, (white triangles) Anomalous mode 
initially excited parallel to probe beam. 



